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Abstract 

We have used the density functional theory to study the effect of molecular elongation on 
the isotropic-nematic, isotropic-smectic A and nematic-smectic A phase transitions of a fluid of 
molecules interacting via the Gay-Berne intermolecular potential. We have considered a range of 
length-to-width parameter 3.0 < xq < 4.0 in steps of 0.2 at different densities and temperatures. 
Pair correlation functions needed as input information in density functional theory are calculated 
using the Percus-Yevick integral equation theory. Within the small range of elongation, the phase 
diagram shows significant changes. The fluid at low temperature is found to freeze directly from 
isotropic to smectic A phase for all the values of xq considered by us on increasing the density while 
nematic phase stabilizes in between isotropic and smectic A phases only at high temperatures and 
densities. Both isotropic-nematic and nematic-smectic A transition density and pressure are found 
to decrease as we increase xq. The phase diagram obtained is compared with computer simulation 
result of the same model potential and is found to be in good qualitative agreement. 

PACS numbers: 64.70.Md, 61.30.Cz, 61.30.Dk 
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I. INTRODUCTION 



Liquid crystal phases which are formed by highly anisotropic complex organic molecules 
have symmetries intermediate between those of isotropic liquid and crystals. The simplest of 
the liquid crystal phases that find application in many electro-optic devices are nematic (N) 
and smectic A (Sm A) ones. In the nematic phase, molecules tend to align along a preferred 
direction called director, breaking the rotational invariance of isotropic liquid (I) but not the 
translational invariance. Partial breakdown of translational invariance along with breaking 
of rotational invariance introduces smectic phases where molecules are essentially confined 
in layers. In Sm A phase molecules are aligned perpendicular to layers with no intralayer or 
interlayer correlation in positions of the center of mass of the molecules. The properties and 
relative stability of these phases are extremely sensitive to the details of molecular structure 
and the true nature of intermolecular interaction potential [1]. Therefore, it is of interest for 
either computer simulation or theory to study the effect of molecular shape anisotropy and 
intermolecular potential on the liquid crystalline phase behaviour and properties. 

Because of the complex structure of the mesogenic molecules it is very difficult to know 
the exact nature of the interaction potential as a function of intermolecular separation and 
orientation. Therefore, modeling of the intermolecular potential only with more physically 
relevant features become inevitable. The pair interaction potential model proposed by Gay 
and Berne (GB)[2] is one such model that includes anisotropic attractive interactions along 
with short range repulsive interactions. It has now become a standard model to study 
liquid crystalline phases and has been widely used in the investigation of various phenomena 
through computer simulation[3-13] and also theoretically [14-18]. 

In the Gay-Berne (GB) pair potential model, the molecules are considered as ellipsoids 
of revolution about the principal axis of the molecule. The interaction potential between 
two such ellipsoidal molecules i and j depends on direction and on the magnitude of 
center-center vector rij = ri — rj and upon molecular axis vectors §i and ej. 
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The GB potential is expressed as 
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The angle dependent range parameters a and strength functions e are given by 
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where do is the smallest molecular diameter and eo is the energy scaling parameter equal 
to the well depth for the cross configuration (ei.ej = r^.ei = r^.e^ = 0). The parameter x 



is a measure of shape anisotropy defined as x 



(Xn 



l)/( x o + 1) w ith iCo being the length- 



to-width ratio of the molecule. Though xo measures the anisotropy of the repulsive core, it 
also determines the difference in the depth of the attractive well between the side-by-side 
and the cross configurations. The parameter x' determines the energy anisotropy defined as 
X 1 = (k' 1 ^ — l)/(k n ^ + 1), where k' is the well-depth ratio for side-by-side and end-to-end 
configuration. The powers \x and v entering in the energy functions are adjustable exponents 
determining the strength of interaction. 

The GB model contains four parameters (xo, k', /i, v) that determine the anisotropy in the 
repulsive and attractive forces in addition to two parameters (o" , eo) that scale the distance 
and energy, respectively. The choices for the values of these parameters are in no way 
unique and they can be varied to yield a wide range of anisotropic potential[3]. The most 



commonly used values of these parameters in the literature are xo = 3.0, k' = 5, fi = 2,v = 1. 
In their molecular dynamic simulation study for this set of parameters de Miguel eta/ [4, 5] 
found three phases, namely, isotropic, nematic and smectic B. No Sm A phase was found to 
be stable. However, existence of Sm A phase has been reported for the set of parameters 
x = 3.0, k' — 5, /i — 1, v — 2 [6]. Recently the effect that changes in the well depth parameter 
k'[7, 15] and molecular elongation xq[8, 9, 17] have upon the overall phase behaviour of the 
system has been investigated. It was found that with the increase in k! the Sm B phase 
is favored at low density while the nematic phase becomes increasingly stable at lower 
temperature as k' is decreased. The variation of molecular elongation parameter xq has 
been found to have a significant effect. An island of Sm A is found to appear in the phase 
diagram for elongations above Xq = 3.0[9]. The range of Sm A extends to both higher and 
lower temperatures as xo is increased. Also the isotropic-nematic (I-N)transition is seen to 
move to lower density (and pressure) at a given temperature with the increase in x . 

In this paper we use the density functional theory (DFT) to study the effect of the 
variation of the length-to-breadth ratio xq on the freezing transitions and freezing parameters 
for the I-N, I-Sm A and N-Sm A transition for a system of molecule interacting via the Gay- 
Berne pair potential. The value of Xq has been varied from 3.0 to 4.0 in steps of 0.2, keeping 
other parameters fixed at k' — 5,// = 2 and v — 1. Pair correlation functions needed as 
input information in the DFT have been calculated using the Percus-Yevick (PY) integral 
equation theory. The paper is organized as follows: In section II we describe in brief the 
density- functional formalism used to locate the freezing transitions and freezing parameters. 
We discuss our results and compare with the available simulation results in section III and 
finally conclude our discussion in section IV. 



II. FUNDAMENTALS OF THE DENSITY-FUNCTIONAL THEORY OF FREEZ- 
ING 

The density functional theory (DFT) directly links the bulk phase behaviour of a fluid 
with its molecular properties. In DFT, the equilibrium density profile of a non uniform 
anisotropic liquid can be determined as the one that minimizes the grand thermodynamic 
potential W regarded as a functional of the single-particle density function p(r, f2) at point 
r and orientation f2. 
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The grand thermodynamic potential has the general form 

- W = (3A - (3p c J rfxp(x) (2.1) 

where A is the Helmholtz free energy, p c the chemical potential and p(x) is a singlet distri- 
bution function, to locate the transition. 

The above mentioned minimum property of the grand thermodynamic potential follows 
from the variational inequality 

W[ Peq ] < W[p] (2.2) 

which is valid for a fixed temperature and chemical potential, where p eq is the equilibrium 
density profile. 

It is convenient to subtract the isotropic fluid thermodynamic potential from W and write 
it as [19] 

AW = W-W f = AW l + AW 2 (2.3) 



with 
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Here Ap(r, fl) = p(r, fl) — pf, where pf is the density of the coexisting liquid. 

The order parameter equation is obtained by minimizing AW with respect to the arbi- 
trary variation in the ordered phase density subject to a constraint that corresponds to some 
specific feature of the ordered phase. This leads to 

P(ri,fli) 

Pf 

where is Lagrange multiplier. Equation(lO) is solved by expanding the singlet distribution 
p(r, fl) in terms of the order parameters that characterize the ordered structures using the 
Fourier series and Wigner rotation matrices . Thus 
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where the expansion coefficients 

21 + 1 r r 
Qi mn (G q ) = -±- J drj dtlp(r,n) x 

exp(-iG q .r)£>^(n) (2.8) 

are the order parameters, G q the reciprocal lattice vectors, po the mean number density 
of the ordered phase and D* n (Q) the generalized spherical harmonics or Wigner rotation 
matrices. Note that for a uniaxial system consisting of cylindrically symmetric molecules 
m = n = and, therefore, one has 

p(r,n) = po££Qi,exp(iG,.r)P,(cos0) (2.9) 
i q 

and 

Qi q = J drj dVtpiv, fi) exp(-iG q .r)P,(cos 9) (2.10) 

where P;(cos#) is the Legendre polynomial of degree / and 9 is the angle between the 
cylindrical axis of a molecule and the director. 

In the present calculation we consider two orientational order parameters 

Pi = Tfffj = (Pi(cos6)) (2.11) 

with 1—2 and 4, one order parameter corresponding to positional order along Z axis, 

fi = Qoo(G z ) = (cos(-z)} (2.12) 

(d, being the layer spacing) and one mixed order parameter that measures the coupling 
between the positional and orientational ordering and is defined as, 

1 2ir 
r = -Q 20 (Gy = (cosi—z^cosO)) (2.13) 
o a 

The angular bracket in above equations indicate the ensemble average. 

The following order parameter equations are obtained by using Eqs.(13) and (17) [18]; 

I rd rn 

Pi — / dz-i I sin9id9iPi{cos9-i) expfs-um] (2.14) 
2d Jo Jo 

p, — — / dz\ \ sin9id9icos( — — ) exp[sum] (2-15) 
2d Jo Jo d 

r = — - / dz\ \ sin9 1 d9 1 P 2 (cos9i)cos( — — ) x 
2d Jo Jo d 

exp[sum] (2-16) 
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and change in density at the transition is found from the relation 



Here 
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and 



cue,) = (^v^E^ + i^ + i^x 

Pi^cosO^CgiliLl; 000) x 

Chu (r 12 )ji (G q r 12 )r 2 12 dr 12 (2.19) 

where Q 1 £j(r) are direct pair correlation function harmonics ,C g (liLl; 000) are Clebsch- 
Gordan coefficients and G g = 27r/d. 
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FIG. 1: Variation of grand thermodynamic potential with smectic interlayer spacing for nematic- 
smectic A transition density p* = 0.219 at T* = 1.40 for the GB potential parameters xq = 4.0, k' = 
5, a = 2 and v = 1. 

In order to evaluate the transition parameters such as order parameters, change in density 
etc, equation (15)-(21) were solved self consistently using values of the harmonics c^uir) 
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evaluated at given temperature and density for each value of xq. The calculation is repeated 
with different values of Sm A interlayer spacing, d. By substituting these solutions in 
Equations (7)- (9) we find the grand thermodynamic potential difference between ordered 
and isotropic phases; i.e. 

- ^ = "AP* + ^P*(2 + Ap*)(7 ° + \{PlCl 2 + 

P*C° U ) + /i 2 ^ + + r 2 C\ 2 (2.20) 

where 

C q LL , = (2L+f) 1 / 2 (2L' + f) 1 / 2 x 

Pf^\ 2 -^) 1/2 C 9 (LL'l-,000)x 

/ c LL/ i(r 12 )ji(G q r 12 )rl 2 dr 12 (2.21) 

J 

At given temperature and density the phase with lowest grand potential is taken as the 
stable phase. Phase coexistence occurs at the value of pf that makes — AW/N = for the 
ordered and the liquid phases. The transition from nematic to the SmA is determined by 
comparing the values of — AW/N of these two phases at a given temperature and at different 
densities. The value of the interlayer spacing d, is found by minimizing the grand potential 
with respect to d. In Fig.l we plot the variation of AW with d*(— d/a ). After selecting 
the value of d for a given density and temperature we locate the transition point using the 
procedure outlined above. In the isotropic phase all the four order parameters become zero. 
In the nematic phase the orientational order parameters P 2 and _P 4 become nonzero but the 
other two parameters p, and r remain zero. This is because the nematic phase has no long 
range positional order. In the SmA phase all the four order parameters are nonzero showing 
that the system has both long range orientational and positional order along one direction. 

III. RESULTS AND DISCUSSION 

The pair correlation functions(PCF's) of molecular fluids are the lowest order microscopic 
quantities which contain information about the structure of the fluid as well as have direct 
contact with the underlying intermolecular interactions. Most of the structural informations 
of an ordered phase are contained in single particle distribution p(x), as shown in the previous 
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FIG. 2: Pressure- Temperature phase diagrams for the GB potential with parameters 3.0 < xq < 
4.0, k' = 5, ji = 2 and z/ = 1 using density-functional theory 



section. In the DFT of freezing the single particle distribution p(x) of an ordered phase is 
expressed in terms of the PCF's of the isotropic fluid. In this study, the PCF's of the isotropic 
phase are found by solving Ornstein-Zernike (OZ) equation using Percus-Yevick (PY) closure 
(details are same as in reference [18]) at different values of the reduced temperatures T*, in 
the range from 0.80 to 1.80 and at various densities for each 3.0 < x < 4.0. 

In Figures 2(a-f) we present the phase diagrams in the pressure-temperature plane. For 
all Xq > 3.0, we first find the direct I-Sm A transition below a certain temperature T* which 
increases from T* ~ 0.99 for x = 3.0 to T* ~ 1.39 for x = 4.0. The lowest temperature 
where Sm A phase starts stabilizing is found to rise as we increase x . On further increase in 
T*, nematic phase enters into the phase sequence. Temperature range of nematic stability is 
found to increase as molecular elongation is increased which is also reflected in the ratio of 
N-Sm A and I-N transition temperatures, T^-a/Ti^n- For example, the value of the ratio 
change from 0.96 to 0.92 as we move from xq = 3.6 to xq = 4.0 with pressure, P*, fixed 
at 3.5. This corresponds to a higher orientational order before the Sm A phase stabilizes. 
According to McMillan theory [20], this value of the T n ^a/Ti^n ratio predicts that the N- 
Sm A transition is first order. Figs. 2(a-f) also show that as x is increased, both the I-N 
and N-Sm A transitions move to a lower pressure. 
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FIG. 3: Variation of I-N and N-Sm A transition densities with molecular elongation parameter xq. 
Dotted and solid lines represent, respectively, I-N and N-Sm A transition densities at T* = 1.40. 
Dot-dashed and long dashed lines are the respective lines at T* = 1.50 and open circles and filled 
circles are those at T* = 1.60, respectively. 

In Figures 3, we plot the I-N and N-Sm A transition densities against the molecular 
elongation x . Both the I-N and N-Sm A transition densities are seen to decrease with the 
increase in xo in agreement to reference [9] . It can be seen from this figure that the difference 
between the I-N and N-Sm A transition densities decreases as x is increased. This is more 
apparent in the low temperature curves where the region between the two transition lines 
gradually shrinks. This makes evident that for a given temperature, with the increase in 
molecular elongation, the Sm A phase become more probable at low density and eventually 
takes over nematic resulting in direct isotropic to Sm A transition. This stabilization of 
smectic phase is to be expected since increase in xq results in deeper well depth for parallel 
configuration of molecules making it energetically more favourable. Also the I-N and N-Sm 
A transition density difference is seen to decrease rapidly as the temperature is reduced at 
a fixed xq, showing the effect of the attractive intermolecular interaction which dominates 
at low temperature facilitating the formation of the Sm A Phase. This also indicates the 
possibility of direct I- Sm A transition for a given x below certain temperature where Sm 
A phase takes over nematic. 

The Sm A phase is characterized by the presence of a non zero translational and mixed 
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FIG. 4: Variation of order parameters with xo for N-Sm A transition at T* = 1.40 keeping other 
GB potential parameters fixed at k' = 5, ji = 2 and v = 1. 

order parameters along with the orientational order parameters. Fig. 4 presents the variation 
of the order parameters at N-Sm A transition with respect to Xq at temperature (arbitrarily 
chosen) T* = 1.40. It shows that the order parameters are monotonically increasing with the 
increase in the molecular elongation. For each xo, we found that the values of translational 
and mixed order parameters, which are characteristics of Sm A phase, decrease rapidly 
with the increase in temperature, eventually becoming zero. This shows that the stable 
Sm A phase becomes progressively less ordered with the increasing temperature and finally 
destabilizes with respect to nematic phase and disappears from the phase diagram. This is 
obvious as attractive forces between the molecules are less important at high temperatures 
and purely repulsive intermolecular interactions do not exhibit Sm A. One should note that 
at x = 3.0, we get stable Sm A phase which has not been observed in computer simulation 
experiments. The phase observed in simulation [9] is Sm B which we have not considered in 
the present study. In Table 1 we have compared our results for the I-N transition parameters 
with the simulation results of refs.[12, 13] at temperature T* = 1.00 for x = 3.0. It can be 
seen that the values of the transition parameters found by us are higher than those of the 
simulations but relatively close to the values of ref. [13]. 

The change in density at transitions Ap*(= Pn ^ pl ) is found to increase as x is increased 
at a given temperature. For example, for T* = 1.50, at I-N transition Ap* changes from 




11 



TABLE I: Comparison of isotropic-nematic transition parameters of the GB(3, 5, 2, 1) potential 
at T* = 1.00. Quantities in reduced units are Pressure P* = Pa^/eo, and /i* = fJ, c / e o- 



Theory 


P* 


A/9* 


P2 


Pa 


p* 




DFT 


0.336 


0.017 


0.675 


0.373 


4.41 


14.57 


MC[13] 


0.320 




0.660 


0.290 






MC[12] 


0.307 


0.017 


0.520 




3.63 


12.79 



4.9% to 6.0% as xq is changed from 3.8 to 4.0, whereas the corresponding change at the 
N-Sm A transition is from 3.4% to 4.7%. In the simulation [9], though the value of Ap* at 
the N-Sm A transition increases from ~ 1.6% to ~ 3.1% as x changes from 3.8 to 4.0, the 
trend at I-N transition is reverse where it decreases from ~ 8% to ~ 5.2%. 

15 r 1 1 1 1 1 1 1 
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FIG. 5: Temperature-Density phase diagram for the GB potential with parameters xq = 3.8, k' = 
5,/x = 2 and v = 1. The solid lines indicate the phase boundaries obtained by using the density- 
functional theory while dashed lines are the simulation results of Brown eta/ [9]. 

In Figs. 5 and 6 we present our phase diagram in the temperature-density plane for 
x = 3.8 and 4.0. The solid lines represent the phase boundaries calculated by us while 
dashed lines are the approximate phase boundaries taken from ref. [9] for the sake of 
comparison. Though the quantitative agreement between the phase diagrams shown in 
the Figures 5 and 6 is not so encouraging, however, they agree qualitatively well. This 
may be due to the fact that in our calculation the transition takes place at higher density 
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FIG. 6: Temperature-density phase diagram for the GB potential with parameters xq = 4.0, k' = 
5, fj, = 2 and v = 1. The curves are the same as in figure 5. 



than those found in the simulations. The critical correlations at which the isotropic phase 
become unstable are found to be higher than actual values because PY theory is known 
to underestimate the orientational correlations [21]. This shortcoming of the PY theory is 
found to increase with increasing temperature which explains why in Figs. 5 and 6 the I-N 
transition boundaries are more shifted than the N-Sm A boundaries. 



IV. CONCLUSION 



We have used the Gay-Berne potential model, to study the effect of the variation of 
elongation parameter x on its phase behaviour with the values of all other parameters 
fixed at k' — — 2 and v — 1. The pair correlation functions of the isotropic fluid 
are calculated using Percus-Yevick integral equation theory and have been used in density- 
functional theory to locate the isotropic-nematic, isotropic-smectic A and nematic-smectic 
A transitions. Within the small range of elongation 3.0 < x < 4.0, the phase diagram 
shows significant changes. We have discussed how various phase boundaries, transition 
densities, pressure, change in density at transition, order parameters etc. change with the 
variation in molecular elongation. We found a stable Sm A phase for all x$ > 3.0. Both 
I-N and N-Sm A transitions are found to move to lower density and pressure as x is 
increased. We have compared our results with those of computer simulations and found 
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that the density-functional theory reproduces all the features of the phase diagrams which 
are in good qualitative agreement. To have a better quantitative agreement there is a need 
to evaluate the isotropic pair correlations more accurately than those given by Percus-Yevick 
theory. Also the use of correlations evaluated directly in the nematic phase may improve 
the results. The work in this direction is in progress. 
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